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Abstract 



The presence of ultra-high-energy cosmic rays (UHECR) results in an increase in the degree of 
ionization in the post-recombination Universe, which stimulates the efficiency of the production of 
^ ' H2 molecules and the formation of the first stellar objects. As a result, the onset of the formation 

'sO , of the first stars is shifted to higher redshifts, and the masses of the first stellar systems decrease. 

As a consequence, a sufficient increase in the ionizing radiation providing the reionization of the 
Universe can take place. We discuss possible observational manifestations of these effects and their 
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■ dependence on the parameters of UHECR. 
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1 INTRODUCTION 

One of the scenarious ^ 01 for the formation of ultra-high-energy cosmic rays (UHECR, with 
' energies above the Greisen-Zatsepin-Kuzmin cut-off, E > 10^^ eV, OEl) suggests they are formed 
c/3 . due to the decay of ultra-heavy X particles with masses > 10^^ GeV. In the interaction of UHECR 
^ ' with low-energy cosmic-microwave background (CMB) photons and the subsequent electromagnetic 
cascades, resonance and ionizing photons are emitted, which can increase fractional ionization of the 
^ . matter in the post-recombination Universe at redshifts ~ 10 — 50 by factors of 5-10 compared to the 
^ [ standard recombination regime. This circumstance could qualitatively change the entire subsequent 
evolution of the Universe, since the degree of ionization substantially influences the rate of radiative 
cooling by gas, and hence, the formation of the first stellar objects. Indeed, cooling of the primordial 
gas at low temperatures T < 10^ K is provided by thermal radiation in rotational lines of H2. In turn, 
H2 molecules can form in the primordial gas only via ion-molecule reactions involving H~ and Hg" , in 
which electrons and protons play the role of catalysts |S]. 

Thus, the presence of UHECR in the early Universe can appreciably influence the subsequent stellar 
phase of its evolution. For this reason, observational manifestations connected with characteristic 
features this phase can be used to constrain the parameters of UHECR. In the present paper, we show 
that due to the influence of UHECR stellar evolution in the Universe begins at earlier epochs, and 
discuss possible observational consequences of this prediction. Section 2 describes the ionization and 
molecular kinetics of baryons in dark halos and their thermodynamics in the presence of cosmic rays. 
In Section 3, we present and discuss our results. A summary is given in Section 4. In calculations we 
assumed a A-CDM model for the Universe: (Qq, f^A, ^^m, f^fe, M = (1-0, 0.71, 0.29, 0.047, 0.72) and 
deuterium abundance n[D]/n = 2.6 x 10~^ [S]. 
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2 COSMIC RAYS AND THE THERMO-CHEMICAL EVOLUTION 
OF BARYONS IN DARK HALOS 



Interaction of UHECR with CMB photons produces high-energy particles - photons, electrons, 
positrons, and neutrinos, which transform through electromagnetic cascades into ionizing Lyc and 
resonant Lya photons at the rate (per unit volume) [HI llUj 

^ = e.A^)H(z)n. (1) 
Here, H{z) is the Hubble constant, n is the baryon density, 

e,,, ^ Miii^Mie^-^eio,, (2) 
1 + z 

is the efficiency of production of ionizing and resonant photons J^, where a is the spectral index of 
the photon spectrum produced in the decay of a superheavy particle, Mig = GeV, Mx is 

the mass of the superheavy particle, Qtot is a function determined by the rate of decay of superheavy 
particles at the present epoch and its redshift dependence [Hl^j. For the cosmic rays to produce 
measurable distortions of the CMB radiation, the function Qtot must be of the order of ~ lO^Mig""'^ 
[H]. Following [7j, we assume that the efficiencies of producing resonant and ionizing photons are 
similar: er{z) = ei{z) = e/(l + z). 

In the standard model of the Universe, in the post-recombination period (z ~ 1100) before the 
formation of the first stars (z ~ 30 — 20), there are no sources of ionizing photons. In these conditions, 
fractional ionization x = ?T,[H"'']/n is determined by photo-recombinations, and ionization by additional 
Lyc and Lya photons produced by the cosmic rays 

± = -k^nx^ + -J—H{z){l - x). (3) 

1 -\- z 

Under such conditions, the rate of formation of molecular hydrogen, / = n[H2]/n, is given by the 
expression 

f = kmnil-x-2f)x, (4) 

where 

™ k3 + k4/{l—x)n k^ + kj / {1 — x)n^ 

is the reduced rate of convergence of H into H2 ^21 Sind ki are the rates of intermediate reactions (see 
the table). At low temperatures (T < 200 K), HD molecules can substantially influence the thermal 
state of the gas ^1 El • The kinetics of these molecules is described by the equation 

g = kmfxndc - nx{kmf + ^02)5, (6) 

where g = n[HD]/n is the number density of HD molecules, dc = n[D]/n is the deuterium abundance, 
and ^Di and A;d2 are the rates of formation and destruction of HD [Tfj . 

The growth of initial perturbations in the dark matter leads to the formation of gravitationally 
bound objects - so called dark halos, in whose potential wells baryons come into virial equilibrium, 
and later on cool radiatively and give birth to stars. 

The role of UHECR at the initial star-formation stage is determined by their influence on the 
thermal state of the baryons. In order to understand this role, we consider the evolution of baryons 
in a virialized halo using the simple system of equations (see, for instance, j2()j ) 

R = u, (7) 
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u — SAj 



(9) 



where n = pj [im^^ pd is the density of dark matter, R is the radius of the region occupied by baryons, 
u is the velocity of the boundary of the baryon cloud, Aj are the cooling heating rates due to Compton 
interactions with the CMB photons and radiation in lines of atomic and molecular hydrogen and HD 
molecules (expressions for the cooling and heating rates are given in [U). The initial values of the 
radius and temperature were assumed to be equal to the virial values for a given halo mass. The 
number densities of electrons and of H2 and HD molecules were calculated for the virialization time, 
assuming a simple prescription for the evolution of the density perturbation JO! • 



Figure 1 shows the z dependences of the temperature, density, and Jeans mass of the gas for three halos 
with the masses of 10'^ Mq, IO^Mq, and lO^M© that have reached a virial state at redshift 2: = 20, and 
for three values of the production efficiency of ionizing photons from cosmic rays (e = 0, 0.1, 1). As an 
example, let us consider in more detail the evolution of a lO^M© halo. In the model with cosmic rays, 
the number density of electrons at the virialization time increases by approximately a factor of five. 
As a result, the number density of H2 molecules substantially increases: during the contraction of the 
baryons, it is 4 x lO"'^ at z = 16.5 for the model with e = 0, while for e = 0.1 it is already 3 x 10"'^. As 
a result, in the models with non-zero e, the gas temperature falls faster, and becomes lower for lower 
densities. A fairly high abundance of molecular hydrogen already at the beginning of the contraction, 
together with low temperature, favore binding almost all the deuterium in HD molecules, due to the 
effects of chemical fractionation |18j . When T < 200 K, the main input to the cooling process is 
provided by HD molecules, and the cooling is then so strong that the temperature drops rapidly, and 
attains the CMB temperature, ~ 50 K. HD molecules provide efficient heat exchange between the 
CMB radiation and baryons via absorption of CMB photons and subsequent transfer of the excitation 
energy to the gas in collisional processes 18^. As a result, the baryon contraction becomes isothermal. 
The gas density grows rather rapidly, and starting from the value > 3 x 10'' cm~^, the optical depth 
in the HD lines exceeds unity. The Jeans mass varies as Mj ~ 100(1 + z)3/2n-i/2M0, and reaches 
1 — 2Mq. At higher densities, three-particle collisions become the major mechanism for forming 
H2, and at n ~ cm the abundance of H2 becomes ~ 1. Thus, the central regions of the cloud 

become opaque in H2 lines [22| [2^] . 

Halos with lower masses evolve much slower, but as the efficiency of producing ionizing photons in- 
creases, the abundances of electrons and H2 and HD molecules also increase, the temperature decreases 
substantially, and the halo begins to contract isothermally but at higher temperatures. For example, 
in the model without cosmic rays, the temperature of a W^Mq halo in this regime is T ~ 300 K, 
which can be explained by the less efficient formation of H2 and HD molecules. Thus, the transition 
of a halo to isothermal contraction is determined by its mass. 

Let us consider the dependence of the redshift zt at which the transition to isothermal contraction 
occurs, on the halo-virialization epoch z^ and the efficiency e, using as an example a halo with mass 
M30- and with a perturbation amplitude corresponding to the 3a level (see, for example, (21]). Figure 
2 shows the dependence zt{zv) for e = 0, 0.1, 1. The upper x axis shows the masses of perturbations 
with amplitudes exceeding 3a that result in the formation of halos with mass of Mh = M30- at z^j-, the 
formation of halos with Mh > M^u has low probability. It is seen that for larger e the halos evolve 
more rapidly to the isothermal regime, particularly in the low-mass end. As the halo mass increases, 
the effect of decrease of the time needed to attain the isothermal state weakens, since cooling and 
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contraction to high densities are possible in massive halos due to energy losses in lines of atomic 
hydrogen. 

Thus, the presence of UHECR substantially accelerates the halo evolution. With increasing e, the 
contraction of the gas proceeds much more rapidly: with e = 0, the baryons in a IO'^Mq halo virialized 
at z = 20 converge to isothermal contraction at z ~ 13.5, while for e = 3 this transition occurs at 
z « 17. The influence of cosmic rays on the thermal evolution can play a principal role for lower-mass 
halos. For instance, when e = 0, halos with ~ 3 x 10^ Mq attain an isothermal state only by 
z ~ 7, when the flux of external UV photons produced by the stellar populations of more massive 
halos becomes sufficiently strong for the temperature of the baryons in such low-mass halos to exceed 
the virial temperature, and for star formation to be suppressed |24| I25| . However, when e = 0.1, 
such halos converge to the isothermal regime much earlier - at z ~ 12, when the influence of the 
ionizing radiation from massive halos is negligibly small. Large values of e result not only in a general 
acceleration of the evolution of halos of all masses, but also in a shift of the minimum mass Mmin for 
the most rapidly evolving halos; this mass corresponds to the maxima of the curves for different values 
of e values in Fig. 2. When e ~ 1, the gas in halos with masses 5 x 10^ attains the isothermal 
regime and conditions favorable for star formation already by z ~ 15 (Fig. 2), i.e. long before the 
influence of earlier formed low-mass halos can become important, while larger-mass objects form later, 
so that their radiation cannot affect the evolution of lower-mass halos. Thus, the influence of UHECR 
is critical for low-mass halos: in the presence of UHECR, the masses of halos in which star formation is 
possible decrease and a larger baryon fraction of the Universe becomes involved into stellar evolution. 
Further, their stellar population can substantially change the ionization and thermal state of the gas, 
since in the hierarchical scenario of galaxy formation low-mass objects contain a larger fraction of the 
mass. This circumstance may turn out to be of a fundamental importance for interpretation of the 
Wilkinson Microwave Anisotropy Probe (WMAP) data 9j. 

According to the WMAP measurements of the background polarization, the optical depth of the 
Universe to Thomson scattering is Te = 0.16 |26j . Such a high value corresponds to the beginning 
of the reionization of the Universe at a redshift of z = 17 26 . On the other hand, the spectra of 
distant quasars show that reionization was completed only by z ~ 6 |271 128j . In addition, the flux of 
ionizing UV radiation expected from quasars and young galaxies is insufficient to reionize the Universe, 
even if the mass distribution of stars, i.e., the initial mass function in the early Universe, was skewed 
towards higher masses. Additional sources of ionizing photons, such as the radiation of accreting black 
holes innilSni, radiation from microquasars |35) and the decay of unstable neutrinos [S21 ESI El] > have 
been recently considered to explain the early reionization. In these conditions, the increase of the 
fraction of baryons condensed into stellar phase at the initial star formation stage in the Universe, 

is an important factor in providing the necessary rate of reionization. In the presence of UHECR, 
this fraction certainly grows, due to the increase in the efficiency of radiative cooling of the baryons. 
Figure 3 shows the dependence of 

on e, where tp{M) is the fraction of baryons present in stars within mass M, calculated using the 
formalism [33] [see (fT^ below], Mmin is the lower limit of masses of dark halos, and M^iji{z, e) is the 
minimum mass of the dark halos that are able to cool and form stars. The corresponding z dependence 
of the number of ionizing photons fuvpp per baryon associated with stellar nucleosynthesis is shown 
in Fig. 3. It is seen from this dependence that at the supposed epoch of reionization of the Universe 
(z < 15), the number of ionizing photons in models with e 7^ can increase by half an order of 
magnitude. Thus, the full reionization of the Universe may be determined by the stimulating effect of 
UHECR on characteristics of the initial stage of stellar nucleosynthesis. 

The presence of UHECR in the early Universe may be revealed from the properties of the char- 
acteristic molecular-line emission arising in the stages preceding formation of the first stars. At the 
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prestellar stage, dense and cold gaseous condensations and molecular clouds form in dark halos, with 
substantial energy losses in H2 and HD lines. The total energy released in rotational lines of H2 and 
HD may be ~ lO^'^erg g~^. The total emission from such molecular clouds in galaxies with masses 
10^ — 10^^ Mq can be substantial, and sufficient to be detected by the planned infrared and submil- 
limeter telescopes ALMA^ ASTRO-F^ and SAFIR^. The influence of UHECR will be manifested as 
an increase in the relative specific energy release in H2 and HD lines 

< >= is...,m)dM ■ 

as well as an increase in the redshift at which this emission is detected, • This dependence is shown 
in Fig. 4. An unambiguous relationship between the observable quantities {Efi^) and ZH2 and the 
parameter e is evident. 

Another physical parameter of the initial stage of stellar evolution that will become accessible to 
observations in the nearest future is the supernova rate. One of the key projects of the James Webb 
Space Telescope^ is to observe supernovae at high redshifts, 2: ~ 10 — 20. Since at least one supernova 
may explode in every low-mass halo the number of supernovae at high redshifts should increase if 
the minimum halo mass is decreased due to the influence of cosmic rays. Figure 5 shows the redshift 
dependence of the expected specific (per unit mass) supernova rate in the Universe for various values 
of e, calculated for a standard hierarchical model of galaxy formation within the formalism [W, for 
the spectrum of perturbations with n = 1. The differential mass distribution of the dark halos at a 
given redshift is 

(M,z)dM=\ —4-- — - — ^ — ^exp ^^—4-7 — rdM (12) 

where Di^z') is the perturbation growth factor, be = 1.69 is the density parameter, and a{M) is the 
rms deviation of the perturbations inside mass M. Summing over mass in (|12)) gives the number of 
objects per unit comoving volume at redshift z. Estimates of the number of supernovae in a halo are 
based on currently existing models 

^(^) = ^ 1.2 X iQ-'^n^M^ + ^)3o'^6 yr-\ (13) 
Tiff 

where (1 + 2)30 = (1 + ^)/^^ and Mq = M/W^Mq. A Salpeter initial mass function is assumed, 
which corresponds to one supernova per every 56M0 = = 0.05^6,5, and = 0.08/^,8 is the 

fraction of baryons which can cool and form stars |HS|. The density of dark matter is p = 200pc = 
200[1.88 X 10~2^/i2(l + z)3], and the free-fah time is tff = (47rGp)"^/2. The efficiency of star formation 
was normalized to the Galactic value, = 0.6%. It is readily seen that at high redshifts, {z > 15), 
the number of supernovae grows approximately linearly with e. 



4 CONCLUSIONS 

We conclude that the increase of the fractional ionization in the post-recombination Universe caused 
by ultra-high energy cosmic rays, has a stimulating effect on the initial stages of stellar evolution in the 
Universe: the epoch of the formation of the first stars is shifted to higher redshifts, and the minimum 
mass for systems in which star formation is possible decreases. This should be manifested in several 
different observational effects sensitive to the efficiency of UV photon production by the cosmic rays, 

^ ht t p : / / www .eso.org/projects/alma/ 
^http://www.ir.isas.ac.jp/ASTRO-F 
^http://safir .jpl.nasa.gov 
*http: / /www. stsci.edu /j wst / 
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e. Considered together, these effects should make it possible to infer the parameters of the UHECR 
theory or obtain additional constraints for them. These manifestations include emission in rotational 
lines of H2 and HD irradiated by contracting protostellar condensations, and the supernova rate and 
its dependence on redshift. Morevover, enhancement of the initial stage of star formation by cosmic 
rays can produce additional ionizing photons, which seem currently to be lacking for providing efficient 
reionization of the Universe. 
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Table 1: Chemical-reaction rates. 
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Reaction coefficient k [cm'^s ^] 
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Figure 2: Dependence of the redshift zt for the 
transition to a state of isothermal gas contraction 
on the epoch of halo formation and the effi- 
ciency e (the solid, dotted, and dashed curves cor- 
respond to e = 0, e = 0.1, and e = 1, respectively). 
The plot shows the evolution of a halo with mass 
M^^ corresponding to 3o" perturbation emerged at 
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Figure 1: Evolution of the (a) temperature, (b) 
density, and (c) Jeans mass for three halos with 
masses of 10^ (dotted), 10^ (dashed), and 10^ 
(solid) that have reached a virial state at redshift 
z = 20, and for three efficiencies of the production 
of ionizing photons from cosmic rays, e = 0, 0.1, 3 
(from left to right). The solid curve in the upper 
plot shows the temperature of the CMB. 
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Figure 3: Redshift dependence of the fraction 
of baryons condensed to a stellar phase in the 
early star- formation stage in the Universe, for 
e = 0, 0.1, 1 (solid, dashed, and dotted, respec- 
tively). A standard model for hierarchical cluster- 
ing is assumed. The scale in the right shows the 
mean number of ionizing (UV) photons emitted 
by the stellar population per baryon, fuupp- 




Figure 4: Redshift dependence of the relative 
specific energy released in H2 and HD lines for 
e = 0, 0.1, 1 (solid, dashed, and dotted, respec- 
tively) . A standard model for hierarchical cluster- 
ing is assumed. 
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Figure 5: Redshift dependence of the supernova 
rate for e = 0, 0.1, 1 (solid, dashed, and dotted, 
respectively). A standard model for hierarchical 
clustering is assumed. 
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